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ABSTRACT 

In the last decade, numerous Lyman-alpha emitters (LAEs) have been discovered with 
narrow-band filters at various redshifts. Recently, multi-wavelength observations of 
I LAEs have been performed and revealed that while many LAEs appear to be young 

pi |. and less massive, a noticeable fraction of LAEs possess much older populations of 

stars and larger stellar mass. How these two classes of LAEs are concordant with the 
hierarchical galaxy formation scenario has not been understood clearly so far. In this 
paper, we model LAEs by three-dimensional cosmological simulations of dark halo 
' merger in a ACDM universe. As a result, it is shown that the age of simulated LAEs 

can spread over a wide range from 2 x lO^yr to 9 x lO^yr. Also, we find that there 
are two types of LAEs, in one of which the young half-mass age is comparable to the 
mean age of stellar component, and in the other of which the young half-mass age is 
appreciably shorter than the mean age. We define the former as Type 1 LAEs and 
the latter as Type 2 LAEs. A Type 1 LAE corresponds to early starburst in a young 
galaxy, whereas a Type 2 LAE does to delayed starburst in an evolved galaxy, as a 
. consequence of delayed accretion of a subhalo onto a larger parent halo. Thus, the 

• ' same halo can experience a Type 2 LAE-phase as well as a Type 1 LAE-phase in the 

' merger history. Type 1 LAEs are expected to be younger than 1.5 x lO^yr, less dusty, 

and less massive with stellar mass Mgtar ^ 5 x 10^ M©, while Type 2 LAEs are older 
than 1.5 X lO^yr, even dustier, and as massive as Mgtar ^ 5 x 10® — 3 x 10^" M©. The 
fraction of Type 2s in all LAEs is a function of redshift, which is less than 2 percent 
at z > 4.5, ^30 percent at redshift z = 3.1, and ^70 percent at z = 2. Type 2 LAEs 
k>( \ can be discriminated clearly from Type Is in two color diagram of z'-H vs J-K. We 

. find that the brightness distribution of Lya in Type 2 LAEs is more extended than 

' the main stellar component, in contrast to Type 1 LAEs. This is not only because 

delayed starbursts tend to occur in the outskirts of a parent galaxy, but also because 
Lya photons are effectively absorbed by dust in an evolved galaxy. Hence, the extent 
of Lya emission may be an additional measure to distinguish Type 2 LAEs from 
Type 1 LAEs. The sizes of Type 2 LAEs range from a few tens to a few hundreds 
kpc. At lower redshifts, the number of more extended, older Type 2 LAEs increases. 
Furthermore, it is anticipated that the amplitude of angular correlation function for 
Type 2 LAEs is significantly higher than that for Type 1 LAEs, but comparable to 
that for Lyman break galaxies (LBGs). This implies that LBGs with strong Lya line 
may include Type 2 LAEs. 

Key words: Galaxies - Lya emitters; Galaxies - Formation; Galaxies - Evolution; 
Galaxies - correlation function 
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1 INTRODUCTION 

To explore the early evolutionary phas es of Ralaxies is impor- 
tant t o understand galaxy formation. [Partridge &: Peebles.1 
* E-mail:shimizu@ipmu.jp l|l967l ') predicted that the starbursts in primeval galax- 

t E-mail:umemura@ccs. tsukuba. ac.jp ies emit significant Lya emission through the recombina- 
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tion of ionized hydrogen in interstellar matter. Although 
many surveys attempted to discover such Lya emitting 
galaxies (Lya emitters: hereafter LAEs), but did not 
succeed to f ind t hem for a long time. In late 1990's, 
ICowie fc Hul (|l998h discovered LAEs with narrow-band fil- 
ters for the first time. Currently, numerous LAEs have 
been discovered at high redshifts 3<z<7by8~ 



10 m class telescopes with narrow-band filters (IHu et al. 


1998 


. 119991. I2OO2I: iKodaira et al.1 l2003l: IShimasaku et al. 


2002 


. 20061: Havashino et al. 2004: Ouchi et al. 


2004. 20051: 


TaniKuchi et al.l 20051: Matsuda et al.l 2004. 2005 


). Although 



the number of observed LAEs increases constantly, the na- 
ture of LAEs is still veiled. Recently, surveys of LAEs in 
the various wavelength bands (opt ical, infrared , sub-mm. 



etc) have been performed activelv ffinkelstein et al. 


2007 


Lai et al."2008': 'Matsuda et al."2007': 'Uchimoto et al. 


2008 


_Finkclstcin ct al. 


2009; Tamura ct al. 2003), and have re- 



vealed that while many LAEs appear to be young and less 
massive, a noticeable fraction of LAEs possess much older 
stellar populations and larger stellar mass. We have not well 
understood how such two classes of LAEs are concordant 
with the hierarchical galaxy formation scenario. As for the 
physical origin of Lya emissi on, the cooling radiation from 
a prirnordial collapsing c loud (|Haiman. Spaans. fc QuataertI 



120001: Fardal et aU 200ll). from a galactic wind- driven shell 



i Taniguchi fc Shiovall200C ), or from star- forming clouds in 



a young starburst galax;y ( Mori. Umemura. fc Ferrarall2004l ) 
has been considered. 

Recently, iMori fc Umemural (|2006l ) proposed a galaxy 
evolution scenario from LAEs to LBGs, based on a 
supernova-dominated starburst galaxy model. In this sce- 
nario, LAEs corres pond to an ear l y evo lutionary phase of 
< 3 X lO^yr. Also. IShimizu et all (|2007l ') have constructed 
an analytic model of LAEs in a ACDM universe, and found 
that if LAEs form in relatively low density regions of the uni- 
verse and the duration of starburst is as short as 0.7 x 10* yr, 
the spatial distributions match the weak angular correlation 
function of LAEs observed at z = 3.1. The spectral energy 
distribution (SED) fitting for observed LAEs has shown that 
LAEs mostly have young average age (~ 10*yr) and low 
stell ar mass (10* ^ 10^ Mr), and are less dus t y or dust 
free ifCawiser et al.ll2006l : iFinkelstein et al 1 l2007l : iLai et al.l 
|2008|). These young LAEs are consistent with the p icture 
bv lMori fc Umemural l|2006l ) and'Shi mizu et al.1 l|2007f ). Very 
recently, deep surveys of LAEs allow us to study detailed 
properties of individual LAEs. As a result, it has been re- 
vealed that LAEs have a wide range of age (10® ~ 10^ yr), 
stellar mass (10® 
Ay up to 1.3 ma g 



10'" Ms 



jFinkelstein et al.ll2007l : iLai et aklliooi : 



iQj, and dust extinction with 

. lagJFin' 

IFinkelstein et alll2009l ). LAEs detected in rest-frame opti- 
cal/near infrared (NIR) bands tend to have older age, larger 
stellar mass, and stronger dust extinction than LAEs un- 
detected in those bands. Thus, the picture of purely young 
starburst galaxies are not always reconciled with observed 
LAEs. So far, the physical reason has not been clarified 
for the existence of an old, massive, and dusty population 
of LAEs. The previous study has shown that a starburst- 
dominated galaxy can emit s trong Lya radiation in d ust-free 
or less dusty environments (|Mori fc Umemurall2006l ). How- 
ever, starburst galaxies cannot be always LAEs in dusty en- 
vironments (e.g. ultra- luminous infrared galaxies). Hence, 
what physical state corresponds to an old population of 



LAEs is an issue of great significance. Some authors argue 
that the clumpy distribution s of dusty gas is important for 
the transf er of Lya photons (|Neufeld.lll99ll : iHansen fc Oh.l 
.2006 : Finkelstein et al.ll2009l ). Since Lya photons undergo 
resonant scatterings on the surface of gas clumps, photons 
can easily escape from the clumpy media. Such an effect 
provides the possibility of old, massive and dusty LAEs. 

In this paper, we explore how a young and old popu- 
lation of LAEs are concordant with a hierarchical galaxy 
formation paradigm. For the purpose, we perform tree- 
dimensional cosmological simulations of dark halo merger 
in a ACDM universe, incorporating the prescriptions of 
star formation, spectral evolution, and dust extinction. 
Throughout this paper, we adopt ACDM cosmology with 
the matter density JIm = 0.3, the cosmological constant 
JIa ~ 0.7, the Hubble constant h = 0.7 in units of Hq = 
100 km s~ ^ Mpc~\ the baryo n density Qsh'^ = 0.02, and 
as = 0.92 (|Spergel et al.ll2003l ). 



2 MODEL AND NUMERICAL METHOD 
2.1 Basic Model 

To pursue the star formation history in the hierarchical 
galaxy formation, we simulate the merging history of sub- 
galactic halos (hereafter subhalo) by three-dimensional cos- 
mological Ai'-body simulations. Here, each particle is re- 
garded as a subhalo that is supposed to consist of dark 
matter and baryons. We simulate A'' = 256'^ subhaloes in 
a comovmg volume of (50 Mpc)^. The mass of a subhalo 
is 2.73 X 1O*M0. It is assumed that the star formation 
is triggered, when a subhalo accretes onto a parent halo. 
Then, we trace the stellar evolution separately for indi- 
vidual subhalos using a spectral synthesis code 'PEGASE' 
((Fioc fc Rocca-Volmerangg Il997l ). Moreover, we take into 
account the effect of dust extinction on Lya emission. The 
present approach allows us to analyze the distributions of 
star forming regions in a halo, and also the clustering prop- 
erties of halos. 

There are basically two types of subhalo accretion. One 
is the almost contemporaneous accretion of subhaloes in a 
young small parent halo, and then coeval starbursts take 
place in the halo. The other is the delayed accretion onto an 
evolved massive halo, and then a newly-triggered starburst 
and an old stellar population coexist. Both types have the 
potentiality of becoming LAEs. If they satisfy LAE condi- 
tions (see the detail below) , we call the former Type 1 LAEs 
and the latter Type 2 LAEs. A schematic view of Type 1 
and Type 2 LAEs is shown in Fig. [T] 



2.2 Numerical Method 

We perform a cosmological A'^-body simulation with 
the particle-particle-partic le-mesh (P^M) algorithm 
(|Hocknev fc Eastwood 198ll). Th e numerical scheme is 
based on I Yoshikawa et al.l (|2000l ) . The size of comoving 
simulation box (Lbox) is set to be the same as the size of 
LAE survey at z = 3.1 by iHavashino et al.1 (|2004l ). that 
is, 50 Mpc in linear scale. This allows us to adjust LAE 
conditions by directly comparing with the observation. 
Here, the periodic boundary condition is imposed. We use 
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Figure 1. Basic conception of a Type 1 and Type 2 LAE. A Type 1 LAE is an early phase of coeval starbursts in a young small halo. 
A Type 2 LAE corresponds to delayed starbursts in an evolved massive halo. 



the Plummer softening function for gravitational force, with 
the softening length of egiav ~ I/box/(10A'^^'^^) (~ 20 kpc in 
a comoving scale). 

A parent halo is fou nd using a friends-of-friends algo- 
rithm ijDavis et al.lll985l ) with linking length equal to 0.2 of 
the mean particle separation. In this study, a system with 
> 2.7 X IO^^Mq (say, equal to or more than 100 particles) is 
identified as a parent halo so that the system mass would be 
corresponding to observed LAEs. As previously mentioned, 
each particle (subhalo) in a parent halo has individual age. 
It is assumed that the star formation is triggered, when a 
subhalo accretes onto a parent halo. The star formation oc- 
curs only in the subhalo, and therefore no star formation is 
triggered at the central of host halo. Also, a subhalo which 
underwent the star formation once does not trigger the star 
formation again. Here, the star formation in each subhalo is 
assumed to occur at the rate as 

i/.(t) = /effexp(--^). (1) 

The star formation timescale Ts is set t o match the typi- 
cal lif etime of young LAE s pred icted by iMori fc Umemural 
l|2006l ') and IShimizu et al.l (|2007l '). that is, 1.0 x 10*yr. The 
efficiency /off is determined so that the final fraction of stel- 
lar mass to total baryonic mass is 10 percent. Here, the 
Salpeter initial mass function is assumed. In the present 
model, each particle is supposed to consist of dark matter 
and baryons. As far as the model LAEs concerned, newly 
accreted subhalos are predominantly responsible for Lyman 
alpha emission, and those subhalos emits the bulk of UV 
photons before their dynamical relaxation. Nonetheless, two- 
component (dark matter and baryon) simulations are very 
worth investigating in the future work. 

We identify LAEs under the same conditions as 
iHavashino et al.1 (|2004), that is, I/Lya.obs ^ 1.4 x 10*^erg/s 
and EWicst ^ 20A, where Z/Lya and EWrcat are observed 



Lya luminosity and Lya equivalent width at rest-frame, re- 
spectively. We estimate observed Lya luminosity by 

-^LycVjObs ~ /csc-£/Lyct,int, (2) 

where I/Lya.int is intrinsic Lya luminosity, and /esc is the 
escape fracti on of Lya photons. Wc calcul ate I/Lyg .int using 
"PEGASE" (jPioc fc Rocca-Volmerangeli 19971 '). The escape 
fraction is evaluated in terms of dust extinction as 

/esc = exp (-Tdust), (3) 

where Tdust is the line-of-sight optical depth of dust with 
opacity in proportion to local metallicity. In practice, Tdust 
is dependent on grain size distributions, which are not well 
known in high-z objects. Hence, we normalize the net cross 
section of dust grains so that the number of simulated LAEs 
should match the obse rved number of LAEs at z = 3.1 
(|Havashino et all |2004| ). (In the present analysis, dust is 
treated as a pure absorber of Lya photons, and the scat- 
tering of Lya photons is not solved. The treatment of dust 
extinction can be more sophisticated in the future analysis.) 

3 RESULTS 

3.1 Two Types of Simulated LAEs 

In Fig. (2] the distributions of mass-weighted age are shown 
for simulated LAEs. The ages spread widely from 2 x lO^yr 
to 9 X 10*yr. Interestingly, a significant number of LAEs are 
much older t han ~ lO^yr that is t he life time of young LAEs 
predic ted by iMori fc Umemural (|2006l ) and [Sh imizu et al.l 
(12003). As seen in Fig. [21 the distributions are fairly con- 
tinuous from younger LAEs to older ones. Younger LAEs 
are early coeval starburst galaxies, while older LAEs result 
from delayed starbursts triggered by later subhalo accre- 
tion onto evolved halos. In order to discriminate delayed 
starbursts from coeval young starbursts, we plot the young 
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Figure 2. Age distributions of simulated LAEs. The age is de- 
fined by the mass-weighted age. 
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Figure 3. Young half-mass ages against mass- weighted mean 
ages. The straight line denotes the equality of two age defini- 
tions. A vertical dashed line shows 1.5 X 10*yr. We define LAEs 
younger than 1.5 X lO^yr as Type Is, and older ones as Type 2s. 
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Figure 4. Equivalent width, Lya luminosity, star formation rate, 
and stellar mass against the mass-weighted age for simulated 
LAEs. A vertical dashed line shows 1.5 X lO^yr. 




half-mass ages against the mass-weighted mean ages in Fig. 
[31 where, the young half-mass age is defined as the mass- 
weighted age of the young half subhalos included in a host 
halo. If starbursts are coeval in a halo, young half-mass ages 
should be comparable to mean ages. It is clearly seen in Fig. 
Othat in the part older than ~ 1.5 x lO^yr, coeval starbursts 
disappear and only delayed starbursts appear to be taking 
place. Hence, in this paper, we define LAEs younger than 
1.5 X 10*yr as Type 1 LAEs, and older ones as Type 2 LAEs. 
It is noted that there is not a clear boundary of two classes 
at 1.5 X lO^yr, but the transition is actually gradual in the 
sense that coeval and delayed starbursts are blended around 
~ lO^yr. Nevertheless, as shown below, we find distinctive 



Figure 5. The fraction of Type 2 LAEs at 1 < ^ < 7. The fraction 
is defined by A^Typo 2 LAEs/^^total LAEs, where A^Type 2 LAEs and 
^total LAEs are the number of Type 2 LAEs and all simulated 
LAEs, respectively. 



trends in photometric properties between Type Is and Type 
2s defined here. 

Fig. 13] shows equivalent width of Lya emission, Lya lu- 
minosity, star formation rate, and stellar mass against the 
mass-weighted age for simulated LAEs. Equivalent width 
(EW) decreases with ages for Type 1 LAEs, ranging from 
40A to 200A, although some are at a level of 350-400A. EW 
for Type 2s is randomly distributed in the range of 30A- 
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150A. Lya luminosity is basically in proportion to star for- 
mation rate for Type 1 LAEs, and gradually decreases with 
ages. For Type 2 LAEs, Lya luminosity randomly spread 
in the range of ~ 2 x 10"'^ erg s~^ to ~ 2 x 10*^ erg s~^. In- 
terestingly, in old Type 2 LAEs (> 6 x lO^yr), a high star 
formation rate does not always lead to high Lya luminosity. 
This can be understood by the effect of dust extinction as 
argued below. 

The stellar mass A/star of LAEs is a fairly mono- 
tonic function of age. Type 1 LAEs are less massive with 
Afstar J; 5 X 10* M©, while Type 2s are as massive as 
Mstar ~ 5 X 10* — 3 X 10^" Mq . Recent observations show that 
LAEs detected by rest-frame optical/NI R bands are more 



mass i ve than 10 Mq and pretty older ( Finkelstein et al 



2007 : iLai et all l2008l: iMatsuda et al.ll2007l : lUchimoto et al. 
2008 ; iFinkelstein et al. 20091 ). Such objects may correspond 
to Type 2 LAEs. 

The fraction of Type 2 LAEs in all LAEs is predicted 
to be a function of redshift. Fig. [5] represents the Type 2 
LAE fraction against redshift. Obviously, the Type 2 frac- 
tion increases with decreasing redshift. At z > 4.5, the Type 
2 fraction is less than 2 percent, since massive haloes have 
not grown yet, whereas it is ~30 percent at redshift z = 3.1 
and ~70 percent at z = 2. This trend is con cordant with the 
recent observations bv lNilsson elal] (|2009l ). 



3.2 Spectrophotometric Properties of Simulated 
LAEs 

In Fig. |6l we show spectral energy distributions (SEDs) of 
a typical Type 1 and Type 2 LAE, wh ich are calculated 
by "PEGASE" l|Fioc fc Rocca-Volmerangeill997l 'l. A Type 1 
LAE is composed of young stars, while the host galaxy of a 
Type 2 LAE is dominated by evolved stars, where a distinct 
4000 A break is seen. Accordingly, the colors of two types 
of LAEs are different. 

In Fig. [71 two color diagrams are shown in observed 
frames for 2 = 2 LAEs and z = 3.1 LAEs. The upper two 
panels are z'-H vs J-K, while the lower two panels are z'- 
K vs J-K, where z'-band is ^9000A, J-band Ril2000A, H- 



Figure 10. Escape fractions of Lyo photons as a function of 
mass- weighted ages of LAEs at z = 3.1. 

band Ril7000A, and K-band Ri22000A. A red line in each 
panel indicates the evolutionary track of stars born in a 
single starburst. Therefore, LAEs near this line represent 
almost coeval starbursts. The mass-weighted ages are shown 
by colors. As expected, Type 1 LAEs 1.5 x 10*yr) well 
follow the single starburst track, whereas colors of Type 2s 
are a function of redshift. At lower redshifts. Type 2s deviate 
farther from the single starburst track. Fig. [7] shows that 
Type 1 and Type 2 LAEs are discriminated more clearly in 
the diagram of z'-H vs J-K, compared to that of z'-K vs J-K. 

3.3 Brightness distribution of a Type 2 LAE 

In Fig. [51 the brightness distributions of Lya emission and 
those at observed-frame K band flux (which corresponds 
to rest-frame 6OOOA flux) are shown for a Type 1 LAE at 
z = 3.1. The angular resolution (pixel size) is set to be 1 arc- 
sec. The brightness distributions at observed-frame K band 
flux trace basically stellar mass distribution. As seen in Fig. 
[51 both distributions are compact (~ 10/i~^kpc), and the ex- 
tents are quite similar to each other. In Fig.[^ the brightness 
distributions are shown for a Type 2 LAE at z — 3.1. The 
Lya emission is diffusely distributed over « 100/i~^kpc or 
stronger at the outskirts, while the brightness distributions 
at observed-frame K band flux exhibit a strong contrast and 
are more concentrated. Moreover, Type 2 LAE is composed 
of some clumps. Such clumpy structures can be seen some 
observed LABs in the SSA22 region ( Uchi moto e t al. 20(51; 
I Webb et al.ll2009l ). This result suggests that some Type 2 
LAEs may not be well dynamically relaxed after mergers, 
since the duration of strong Lya emission is shorter than 
the relaxation time of such a large system. 

Since Lya emission is radiated mostly by starbursts, 
Lya bright regions correspond to sites of subhalo accretion. 
These results suggest that if we observe LAEs in various 
wavelengths, the extent and morphology of Type 2 LAEs 
appear to be different in each band. Recently, such an off- 
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Figure 7. Two color diagrams in observed frames. The upper two panels are z'-H vs J-K, where the left panel is shown for LAEs at 
z = 2 and the right panel is at z = 3.1. The lower two panels are z'-K vs J-K, where the left panel is shown for LAEs at z = 2 and the 
right panel is at z = 3.1. A red line in each panel denotes the evolutionary track of a single starburst. The mass- weighted ages are shown 
by colors with an attached color legend bar. 



set has been reported in local starburst galaxy (|Haves et al.l 
l2007l : lOstlin et al.ll2009ll . whil e no offset is found in LAEs 
at a higher redshift of 2 = 4.4 l|Finkelstein et al.ll20oi ). This 
can be related to the result that the number of more ex- 
tended, older Type 2 LAEs increases at lower redshifts, as 
shown below. 

Lya emissions are affected more severely by dust extinc- 
tion than continuum radiation. The effect of dust extinction 
is expected to be stronger for old metal-enriched systems. 
In Fig. IIOI the escape fractions of Lya photons are plotted 
against ages of LAEs. The escape fractions are a decreasing 
function of LAE ages as anticipated. In particular, in Type 
2 LAEs, the escape fractions decrease down to « 10 percent. 
Thus, for aged Type 2 LAEs, Lya equivalent widths or lu- 
minosities become smaller than expected purely from star 
formation rate (see Fig.|4]). Nevertheless, in the outskirts of 
halo, starbursts by delayed subhalo accretion take place fre- 
quently, where the dust extinction is relatively small. Hence, 
Lya emissions can be brighter in the outskirts of halo. 

As shown in Fig. [9l the size of Type 2 LAEs can ex- 
ceed 100 kpc in physic al size. This is compa rable to the size 
of Lya blobs (LABs) (jlVIatsuda et al.ll2004l ). Thus, some of 
Type 2 LAEs with EW>100A may account for a part of 
observed LABs. 



3.4 The halo-size distributions of simulated LAEs 

Here, we analyze the halo-size distributions of simulated 
LAEs. Fig. [11] represents the halo-size distributions of all 
simulated LAEs as a function of mass wighted age at z = 
2,3.1, and 4.5. Here, the radius from the center of gravity 
of a halo within which 95% of the total mass is included is 
defined as the size of a halo. At z = 4.5, there are no Type 
2 LAEs. At lower redshifts, Type 2 LAEs appear and the 
number of more extended, older Type 2 LAEs increases with 
decreasing redshifts. Interestingly, the range of Type 2 LAE 
size is quite broad. Small Type 2 LAEs with the size of a few 
tens kpc are as compact as Type 1 LAEs, while large Type 
2 LAEs with the size of a few hundred s kpc are comparable 
to Lya blobs (LABs) ijMatsuda et al.l[2004ll . 

3.5 Clustering Properties of Two-type LAEs 

In order to explore the clustering properties of two types of 
LAEs, we calculate a two-point angular correlation function 
(ACF) and a two-point angular cross-correlation function 
(CCF) between Type 1 and Type 2 LAEs. Fig. [T2lrepresents 
ACFs of each type LAEs and all simulated LAEs. Also, the 
ACF of LAEs observed in the SSA22 field l|Havashino et all 
l2004h is shown. ACF of Type 1 LAEs as well as ACF of all 
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Figure 8. The brightness distributions of a Type 1 LAE at z = 3.1. Left and right panels show the brightness distributions of hya 
emission, and those at observed-frame K band flux (which corresponds to rest-frame 6000A fiux), of stellar components, respectively. 
Each color bar shows the flux of hya emission and observed-frame K band flux, respectively. The angular resolution (pixel size) is set to 
be 1 arcsec. 
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Figure 9. Same as Fig.|8l but for a Type 2 LAE at z = 3.1. 
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Figure 11. The halo-size distributions of simulated LAEs as a 
function of mass weighted age at z = 2, 3.1, and 4.5. 
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Figure 13. Two-point angular cross-correlation function (CCF) 
of simulated LAEs. CCF between Type 1 and Type 2 LAEs is 
depicted by a dash-dot line. For comparison, ACF of all simulated 
LAEs is also shown by a solid line. 



detected as SMGs, our simulation is consistent with this ob- 
servation. 
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Figure 12. Two-point angular correlation function (ACF) of all 
simulated LAEs and each type LAEs at z = 3.1. A solid line is 
ACF of all simulated LAEs, and dashed and doted lines show 
ACF of Type 1 and Type 2 LAEs, resp ectively. Open-squares are 
the ACF of LAEs observed in SSA22a l lHavashino et al.ll20o4) . 



simulated LAEs is quite weak and well matches the observed 
ACF in in the SSA22 field, whereas ACF of Type 2 LAEs is 
significantly stronger than that of Type 1 LAEs. This implies 
that Type 1 LAEs prefe rentially reside lower- density regions 
at z = 3.1, as shown bv lShimizu et all |20o3), while Type 2 
LAEs are located in higher-density regions. The correlation 
strength of Ty pe 2 LAEs is comparab le to that of LBGs in 
SSA22 region (|Giavalisco et al]fl998l ). This result suggests 
that Type 2 LAEs may be a subsample of LBGs with strong 
Lya emission. 

Fig. [THl shows CCF between Type 1 and Type 2 LAEs. 
The CCF is st ill weaker than ACF of aU simulated LAEs. 
Very recently, iTamura et al.l l|2009l) have found that the 
cross correlation between Sub-mm galaxies (SMGs) and 
LAEs in the SSA22 field is very weak. Type 2 LAEs pos- 
sess a large amount of dust. If many Type 2 LAEs can be 



4 CONCLUSIONS AND DISCUSSION 

To explore the origin of two populations of LAEs recently 
found, we have performed three-dimensional cosmological 
Ai'-body simulations of subhalo merging history in a ACDM 
universe. We have incorporated star formation history, SED 
evolution, and dust extinction. As a result, we have found 
that the age of simulated LAEs can spread over a wide range 
from 2 X lO'^yr to 9 x 10**yr. Also, we have revealed that 
there are two types of LAEs. We have defined LAEs younger 
than 1.5 x lO^yr as Type Is, and older ones as Type 2s. In 
Type 1 LAEs early coeval starbursts occur in small parent 
halos, while in Type 2 LAEs delayed starbursts take place in 
evolved massive haloes as a consequence of delayed accretion 
of subhalos. A parent halo can experience repeatedly a Type 

2 LAE-phase after a Type 1 LAE-phase. 

The stellar mass of Type 1 LAEs is Mstar < 5 x 10** M©, 
while Type 2 LAEs are as massive as Mstar ~ 5 x 10* — 

3 X 10^" Mq. The physical properties of Type 1 and Type 
2 LAEs are concordant with those of t wo populations of 



LAEs observed with multi-wavelengths ( 


Finkelstein et al. 


I2OO7!: iLai et al.l 12008"; 'Matsuda et al.l 12007 


: Uchimoto et al. 



2008; Finkelstein et al. 2009). The fraction of Type 2s in aU 
LAEs is a function of redshift, which is less than 2 percent at 
z > 4.5, ~30 percent at redshift z = 3.1, and ~70 percent at 
z = 2. Th is trend is con si stent with two populations of LAEs 
found by iNilsson et al.1 (120091 '). Type 2 LAEs are expected 
to be discriminated clearly from Type 1 LAEs in two color 
diagram of z'-H vs J-K. We find that the brightness distri- 
bution of Lya in Type 2 LAEs is more extended than the 
main stellar component, in contrast to Type 1 LAEs. This is 
not only because delayed starbursts tend to occur in the out- 
skirts of a parent galaxy, but also because Lya photons are 
effectively absorbed by dust in an evolved galaxy. The sizes 
of Type 2 LAEs range from a few tens to a few hundreds 
kpc. At lower redshifts, the number of more extended, older 
Type 2 LAEs increases. Small Type 2 LAEs are as compact 



Two types of LAEs 9 



as Type 1 LAEs, while large Type 2 L AEs exceediriR 100 kp c 
are comparable to Lya blobs (LABs) (jMatsuda et al.ll2004 ). 

Moreover, we have found that the clustering of Type 2 
LAEs are even stronger than Type 1 LAEs. The amplitude of 
angular correlation function of Type 2 LA Es is comparable 
to th at of Lyman break galaxies (LBGs) (|Giavalisco et al] 
1 19981 ). This suggests that LBGs with strong Lyo? line can 
be Type 2 LAEs. The two-point angular cross-correlation 
function is still weaker than that of all LAEs. If many Type 
2 LAEs can be detected a s SMGs, this result i s consistent 
with recent observation bv lTamura et al.l (|2009l ). 

Interestingly, in a low redshift universe at 0.2 < z < 
0.35, the Galaxy Evolution Explorer (GALEX) have found 
many LAEs older than 2 x 10*yr, more massive than 10^ 
in stellar compo nent, and having smal l escape fractions 
of Ly g photons jPeharveng et all l2008l : iFinkelstein et all 
I2OO9I ). These properties are quite similar to those of Type 2 
LAEs we defined in this paper. Hence, local LAEs could be 
a good sample to study the detailed physical states of Type 
2 LAEs. 

Furthermore, in an evolved population of LAEs, active 
galactic nucleus (AGN) events may be anticipated, since su- 
permassive black holes reside evolv ed galaxies in a local uni- 
verse fe.g .. iMarconi fc Huntll2003l . and references therein). 
Recently, lYamada et al. ( 20091 ) have found that AGNs are 
exclusively (90%) associated with the massive objects with 
stellar mass larger than lO^'''^ Mq. Also, GALEX ioMnA that 
in lo w-redshift LAEs an AGN fraction is as high as forty per- 
cent (|Finkelstein et al]l2009l ). These findings may imply that 
Type 2 LAEs could have a higher fraction of AGNs. Hence, it 
seems worth investigating the contribution of AGNs to Lya 
emission especially for Type 2 LAEs (Shimizu & Umemura, 
in preparation). 
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